Prediction of polar ordered oxynitride perovskites 
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Abstract 

Using a materials by design approach, we find a class of ordered oxynitride piezoelectrics with 

d 

perovskite structure. We predict that ordered YSiC^N and YGeC^N are characterized by large 
nonlinear optic responses and by some of the largest polarizations known to-date. 
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Most piezoelectrics (PZs) or polar dielectrics are oxides, and those commonly used are 
Pb-bearing compounds. PZs transform mechanical energy into electrical energy and vice 
versa, and serve many uses ranging from medical ultrasound and actuators to non-linear 
optic devices [Tj. Materials for high temperature and/or high electric field applications are 
highly desirable, like Pb-free PZ materials as mandated by European law. 

We have applied a synthetic solid state chemistry approach using ab initio density func- 
tional calculations [21 El HI El O [7] to search for useful materials. Instead of an automated 
search of hundreds of compositions, we used a chemical approach similar to what would 
be done to experimentally search for materials. We built the oxynitride structures starting 
from the ideal cubic ABO3 perovskite structure by replacing in each unit cell one oxygen 
atom by one nitrogen (Fig. 1). In this way we aim to drive a stronger polar distortion due 
to the strong covalent bonding expected between the octahedrally coordinated B-cation and 
N, as well as the higher nominal ionic charge of N 3_ compared with O 2- . The resulting 
structures are ordered and polar with tetragonal symmetry, PAmm space group and five 
atoms per unit cell. In order to balance the -7 valence of O2N we tried several combinations 
of trivalent and tetravalent cations on the respectively A and B sites: YSi0 2 N, YGe0 2 N, 
InTi0 2 N, GaTi0 2 N, BiTi0 2 N, YTi0 2 N, BiZr0 2 N, YC0 2 N, YSn0 2 N and YZr0 2 N. 

We used the local density approximation (LDA) of the density functional theory [H |9] with 
a planewave basis and pseudopotentials in the ABINIT implementation [10]. We employed 
Troullier-Martins pseudopotentials [UJ with a 40 Ha kinetic energy cut-off. We computed 
the dynamical and dielectric properties and the coupled response of the energy to the stresses 
and electric fields using density functional perturbation theory [121 E2 EH]- 

We relaxed the PAmm atomic positions and cell parameters at zero pressure for each 
of the above compositions and checked the insulating/metallic character computing the 
electron band structure and the corresponding density of states. We found that InTi0 2 N 
and GaTi0 2 N are metallic and all the other compounds are insulators. We checked the 
dynamical stability for each composition computing the phonon band dispersion and from 
all the combinations of cations listed above only YSi0 2 N and YGe0 2 N were dynamically 
stable. 

The crystal structure of the PAmm modification of YSi0 2 N and YGe0 2 N is shown in 
Figure [T] and the relevant structural parameters are listed in Table |T] All the atoms are 
displaced relative to their positions in the centrosymmetric structures. The corresponding 
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cation displacements are 0.76 A for Y and 0.61 A for Si in YSi0 2 N and 0.82 A for Y and 
0.77 A for Ge in YGe0 2 N. These large displacements enhance the polar character of the 
structures but also deepen the potential wells: 1.54 eV per unit cell in YGe0 2 N and 0.63 
eV in YSi02N. Due to these large potential wells, the polar structures will probably not be 
switchable, thus these are not predicted to be ferroelectric without modification. 

Oxynitrides are generally produced by low-temperature ammonization processes, which 
do not allow chemical ordering to occur. [151 HU [TTJ [18]. We have tested different atomic 
configurations and found that the ordered P4mm structure is the lowest in energy, e.g. lower 
by at least 0.25 eV/0 2 N group for YSi0 2 N than any other disordered configuration. We 
conclude that even if ordering in the oxynitrides is difficult [32] , the previous experimental 
structures were disordered because of fabrication techniques. We suggest that the ordered 
materials should be producible by high temperature methods or by molecular beam epitaxy. 
Thus here we discuss ordered polar materials. 

The phonon band dispersions for the two structures are shown in Figure [2] All the modes 
are stable, indicating dynamical stability and thermodynamical (meta) stability. The elastic 
constants also indicate mechanical stability (Table The heat of formation of the oxyni- 
trides relative to a mixture of YN (fee structure) and Si02 or Ge02 (quartz structure) shows 
stability for the PZ phases (-0.387 eV/unit cell for YSi0 2 N and -0.309 for YGe0 2 N). With 
these stability criteria the PZ structures of the two oxynitrides are at least thermodynam- 
ically metastable and their large negative heat of formation suggests that they should be 
easily synthesized. 

Both materials are insulators with large LDA band gaps, 2.6 eV in YGe0 2 N and 4.4 eV 
in YSi0 2 N. The computed dielectric tensors are e°°=[4.37 4.37 4.49] and e°= [13.44 13.44 
9.71] for YSi0 2 N and e°° = [4.60 4.60 5.01] and e°=[16.10 16.10 9.80] for YGe0 2 N. Because of 
symmetry, the Born effective charges, Zf* = 1/Q SPi/Sr", where Q is the unit cell volume, 
Pi is the polarization along the direction i and r° is the displacement of atom a along 
direction j at zero electric field, are diagonal. Their values in the PZ structures are listed in 



Table III The Born effective charges are only slightly enhanced due to the absence of the 
p-d hybridization characteristic to other oxide perovskites [201 \21\ . 

We computed the spontaneous polarization using the modern theory of polarization 
[22] [23l 124"] . The Berry's phase polarization is a lattice of values, with a lattice spacing 
of 2eR/fi, where e is the electron charge and R are the lattice vectors. Polarization differ- 
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ences, however, are unambiguous. Experiments measure polarization changes (usually via 
hysteresis loops) and absolute polarizations are not measurable using standard methods. For 
ferroelectrics with small polarizations comparison between theory and experiment is straight- 
forward, but in materials with large polarizations the lattice of polarizations measured for 
a single polar structure cannot be unambiguously reduced to an experimentally meaningful 
value, called the effective polarization, without additional computations [26] . For tetrag- 
onal symmetry, only the z-component of P is non-zero. We obtained the formal polarization 
along the polar clXIS Z cLS P Bevvy — 

130(±n308)/iCy cm 2 for YSi0 2 N and 103(±n293)/xC/cm 2 
for YGe0 2 N, where the polarization lattice spacing is indicated in parentheses. In order to 
find the effective polarization, we computed its value using the Born effective charges, Z* (Ta- 



ble III) and the displacement vectors from the ideal perovskite structure, u: AP = YJ Z*u. 
This gives effective polarizations Pz=-163 /iC/cm 2 for YSi0 2 N and P^=-171 /zC/cm 2 for 
YGe0 2 N. These values are not identical to the Berry's phase because of the linear approx- 
imation, but are very close. Using these values, we see that the effective P Berry is -178 
/xC/cm 2 for YSi0 2 N and -190 /zC/cm 2 for YGe0 2 N and these are our best estimates. The 
direction of polarization is from the B ion towards the closer N. Figure [3] shows the polar- 
ization lattices and comparison of the effective polarization and the Berry's phase values. 
The Born effective charges vary with distortion (Fig. [3^ and d), the absolute values of the 
charges being slightly larger for the centrosymmetric phase than the ferroelectric phase. 

These values of the polarization are amongst the highest ever reported so far in the 
literature. For comparison the spontaneous polarization is 71 /iC/cm 2 in LiNbO"3, 50 in 
LiTa0 3 ,[27j, 70 in Pb(Mg,Nb)0 3 [2E] and >150 for polycrystalline thin films of BiFe0 3 [29]. 
Other large polarization are those computed for pure PbTi0 3 (88) |3Q] . BiFe0 3 (90-100) [25J, 
PbV0 3 (150) [31] • 

We computed the piezoelectric constants tensors, d^ = de^jde^ where e is the stress 
tensor and e the electric field, using the linear response theory of elasticity [HI [32] and 
obtained d m =-12.5 pC/N, d 333 =-8.3 pC/N and dn 3 =0.4 p/N values for YSi0 2 N and d m =- 
20.5 pC/N, d 333 =-5.5 pC/N and d 113 =-0.5 pC/N values for YGe0 2 N. We also calculated 

(2) (2) 

the non-linear optical coefficients, dijk = XijL/2 where xljk = 3E £i£: > £k /V is the third-order 
derivative of the energy with respect to electric fields, via linear response [33], and ob- 
tained dn 3 =2.03 pm/V and d 333 =-5.5 pm/V in YSi0 2 N and dn 3 =2.63 pm/V and d 333 =- 
4.51 pm/V in YGe0 2 N. The highest value to-date is shown by AANP, an organic material 
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with d3n=dn3=80 [34J. In KH2PO4, a commonly used non-linear optical material, the non- 
linear optical coefficients are on the order of 0.4 pm/V and in ADP on the order of 0.8 pm/V 
[55] . The electro-optic coefficients, c^, obtained from A(e -1 )y = J2k( c ijk^k) where (e _1 )r/ 
is the inverse of the dielectric tensor, have values of Ci33=-1.43 pm/V, C333 =0.64 pm/V and 
c 33i=-l-32 pm/V for YSiC^N and ci 33 =-0.84 pm/V, 0333=1. 06 pm/V and C33i=-1.77 pm/V 
for YGe0 2 N. 

The high spontaneous polarizations, the non-linear optic coefficients, and the large stabil- 
ity of the piezoelectric structures make both YSiC^N and YGeC^N promising materials with 
broad possible applications as piezoelectrics and non-linear dielectrics. In addition, they may 
be useful as X-ray and neutron generators via accelerations of ions at their surfaces [36J due 
to their high polarizations. 

Acknowledgments 

Razvan Caracas acknowledges I. Grinberg for useful discussions. We acknowledge support 
from the U.S. Office of Naval Research. 



5 



[1] A. Arnau (Ed.), Piezoelectric transducers and applications, Springer Verlag, Heidelber (2004). 
[2] A. Franceschetti and A. Zunger, Nature 402, 60 (1999). 
[3] S.C. Erwin and I. Zutic, Nat. Mater. 3, 410 (2004). 

[4] P. Baettig, C.F. Schelle, R. LeSar, U.V. Waghmare and N.A. Spaldin, Chem. Mater. 17, 1376 
(2005). 

[5] S.V. Dudiy and A. Zunger, Phys. Rev. Lett. 97, 046401, (2006). 
[6] D.I. Bile and D.J. Singh, Phys. Rev. Lett. 96, 147602 (2006). 
[7] K. Kang, Y.S. Meng, J. Breger, CP. Grey and G. Ceder, Science 311, 977 (2006). 
[8] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 
[9] W. Kohn and L.J. Sham, Phys. Rev. 140, A1133 (1965). 
[10] X. Gonze, J.-M. Beuken, R. Caracas, F. Detraux, M. Fuchs, G.-M. Rignanese, L. Sindic, M. 
Verstraete, G. Zerah, F. Jollet, M. Torrent, A. Roy, M. Mikami, Ph. Ghosez, J.-Y. Raty and 
D.C. Allan, Comput. Mater. Sci. 25 478 (2002). 
[11] N. Troullier and J.L. Martins, Phys. Rev. B 43, 1993 (1991). 

[12] S. Baroni, S. de Gironcoli, A. Dal Corso and P. Giannozzi, Rev. Mod. Phys. 73, 515 (2001). 
[13] X. Gonze, G.-M. Rignanese and R. Caracas, Z. Kristall. 220, 458 (2005). 
[14] D.R. Hamman, K.M. Rabe and D. Vanderbilt, Phys. Rev. B 72, 033102 (2005). 
[15] S.J. Clarke, B.P. Guinot, C.W. Michie, M.J.C. Calmont and M.J. Rosseinsky, Chem. Mater. 
14, 288 (2002). 

[16] Y.-I. Kim, P.M. Wodward, K.Z. Baba-Kishi and C.W. Tai, Chem. Mater. 16, 1267 (2004). 
[17] R. Marchand, Y. Laurent, J. Guyader, P. L'Haridon and P. Verdier, J. Eur. Ceram. Soc. 8, 
197 (1991). 

[18] M.S. Hegde, Proc. Indian Acad. Sci. (Chem. Sci.) 113, 445 (2001). 

[19] B. Ravel, K.-I. Kim, P.M. Woodward and CM. Fang, Phys. Rev. B 73, 184121 (2006). 

[20] R.E. Cohen, Nature 358, 136 (1992). 

[21] Posternak, M., R. Resta, and A. Baldereschi, Phys. Rev. B 50, 8911 (1994). 
[22] R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 1651 (1993). 
[23] D. Vanderbilt and R. D. King-Smith, Phys. Rev. B 48, 4442 (1994). 
[24] R. Resta, Rev. Mod. Phys 66, 899 (1995). 



6 



[25] J. B. Neaton, C. Ederer, U.V. Waghmare, N.A. Spaldin, and K.M. Rabe, Phys. Rev. B 71, 
014113 (2005). 

[26] R. Resta and D. Vanderbilt, in Physics of Ferroelectrics: a Modern Perspective, Eds. C.H. 

Ahn, K.M. Rabe, and J.M. Triscone, Springer- Verlag, New York (2007). 
[27] M.E. Lines and A. M. Glass, Principles and applications of ferroelectrics and related materials, 

Oxford Classics, Oxford University Press, Oxford (2001). 
[28] Z. Kutnjak, J. Petzelt and R. Blinc, Nature 441, 956 (2006). 

[29] K.Y. Yun, D. Ricinschi, T. Kanashima, M. Noda and M. Okuyama, Japan. J. Appl. Phys. 43, 
L647 (2004). 

[30] G. Saghi-Szabo, R.E. Cohen and H. Krakauer, Phys. Rev. Lett. 80, 4321 (1998). 
[31] U. Yoshitaka, S. Tatsuya, I. Fumiyuki and O. Tamio, Jpn. J. Appl. Phys. 44, 7130 (2005). 
[32] Z. Wu and R.E. Cohen, Phys. Rev. Lett. 95, 037601 (2005). 
[33] M. Veithen, X. Gonze and Ph. Ghosez, Phys. Rev. B 71, 125107 (2005). 
[34] S. Tomaru, S. Matsumoto, T. Kurihara, H. Suzuki, N. Ooba and T. Kaino, Appl. Phys. Lett. 
58, 2583 (1991). 

[35] R.C. Eckardt, R.L. Byer, H. Masuda and Y.X. Fan., IEEE J.Q. Electr. 26, 922 (1990). 
[36] B. Naranjo, J.K. Gimzewski and S. Putterman, Nature 4343, 1115, (2005). 



7 



TABLE I: Structural parameters for YSi0 2 N and YGe0 2 N in the P4mm ferroic structures. Y, 
Ge/Si, and N atoms are in lb(l/2 1/2 z), la(0 z), 2c(l/2 z) and la Wyckhoff positions 
respectively. 

a(A) c(A) zy z Ge/Si z z N 
YSi0 2 N 3.228 4.435 0.342 0.889 0.027 0.513 
YGe0 2 N 3.307 4.66 0.325 0.874 0.038 0.500 



TABLE II: Elastic constants (in GPa) at zero pressure computed from density functional pertur- 
bation theory [El Q3] . 

11 ^33 W2 W3 W4 ^66 

YSi0 2 N 466 317 211 100 135 231 
YGe0 2 N 374 345 176 86 90 167 



TABLE III: Born effective charge tensors, Z*- in the ferroelectric phases of YGe0 2 N (columns 1-4) 
and YSi0 2 N (columns 5-8). Due to symmetry the tensors are diagonal. 
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FIG. 1: Theoretical P4mm polar structure of the ordered ABO2N oxynitride perovskites. The 
most promising materials we found have Y in the A site, and Si or Ge in the B site. 
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YGe0 2 N 




FIG. 2: Phonon band structure computed in the density functional perturbation theory for YSiO^N 
(a) and YGeO^N (b). The path through the Brillouin zone passes through: T=(0 0), X=(l/2 
0), M=(l/2 1/2 0), Z=(0 1/2), R=(l/2 1/2), A=(l/2 1/2 1/2). 
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FIG. 3: The lattice of polarizations for YGeO^N (a) and YSiO^N (b) and the Born effective 
charges for YGeO^N (c) and YSi02N (d) as a function of the position (=distortion) along the 
centrosymmetric-to-ferroelectric direct path. 
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